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Krstulovic and Brachet reply: In the preceding 
Comment pQ , Kozik raised a criticism against the bottle- 
neck proposed in our Letter [2] to be causing a thermal- 
ization delay when dispersive effects, controlled by the 
coherence length £, are large at truncation wavenumber: 
&max- The late-time energy spectrum presents a front at 
wavenumber k c (t) propagating toward higher wavenum- 
bers and leaving in its wake a quasi-thermalized distri- 
bution. Kozik argues that our observations agree with 
the relaxation scenario, developed by Svistunov [3], that 
involves no bottleneck and predicts k c (t) ~ i 1 / 4 . 

Indeed, it is apparent on FigJT] (where k c (t) ~ t a corre- 
sponds to a line of slope (a—l)/a) that four out of eleven 
runs (vi, vii, viii and xii) are somewhat compatible with 
the Svistunov prediction. However, the prediction only 
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FIG. 1: (color online). Parametric representation dk c /dt v.s. 
fcc/^max (adapted from Fig.4.e of [2]). Different scaling laws 
are displayed. Svistunov prediction: solid line (fc c /fc max )~ 3 ; 
fit to run xi: large dashed line (fc c /fcmax) -2 ' 4 ; dotted line 
(fcc/fcmax) 1 and dotted-dashed line (fc c /fc max ) _ °°. 



works in the limited range 0.4 < k c /k max < 0.8 where run 
xi, with Taylor-Green initial data and £fc max ~ 6, yields 
a slope of —2.4. Runs vi, vii and viii, with initial data 



prepared using the stochastic Ginzburg-Landau equation 
and £/c max ~ 6, have slopes closer to the Svistunov pre- 
diction of —3. Runs that saturate with k c /k ma x ~ 1 
are reaching (truncated) thermal equilibrium and are not 
spectrally well-converged. In contrast, runs i-iv, with 
£fcmax ~ 24, saturate at fc c /^max < 0.4 and are well- 
converged but the data suggests a logarithmic growth of 
k c (t) (vertical line on Fig[T]), a behavior very different 
from that predicted in [3J. 

This discrepancy is perhaps due to the fact that Svis- 
tunov considers a two stages process: first a condensation 
produced by a particle-flux wave propagating to low en- 
ergies and then a wave propagating from the low to high 
energy region. It is not absolutely clear that the initial 
conditions of our Letter [2J, really correspond to any of 
the stages considered by Svistunov (see the discussion 
following Eq. (4.7) of [3]). 

Concerning the criticism against our use of the word 
"bottleneck", we believe it is related to a limitation in 
Svistunov theory. Indeed, it is well known that Bogoli- 
ubov's dispersion relation u>B(k) = fcc(l + fc 2 £ 2 /2) 1//2 
(where c is the sound velocity) implies (around wavenum- 
ber k ~ l/£) a change from propagative to dispersive 
behavior. This elementary point is not completely ad- 
dressed in Svistunov theory, in particular at level of the 
kinetic equations 3.10-3.13 of [3] and Eq.(l) of p] [?]• 
Thus Svitunov's analysis is only applicable for wavenum- 
bers k >> l/£. This limitation does not allow one to 
appreciate the importance of £ and to grasp that fc£ (in 
particular £,k max ) is an important dimensionless param- 
eter in this problem leading to a crossover between dif- 
ferent regimes (see Fig|l]and also Fig.7b-c of [1]). 

In a physical BEC, fc ma x correspond to the equipar- 
tition wavenumber fc Gq (see [2] and Sec. IV of [4]). 
Sinatra and Castin j5] have shown that the slowdown 
of thermalization reported in [5] can be related to the 
behavior of the (classical) damping rate around equilib- 
rium that reaches a maximum around fc£ ~ 3 and de- 
cays for k£ >> 1. They have established that, at fixed 
k^ well beyond its maximum, the (quantum) Beliaev- 
Landau damping rate approaches the classical one pro- 
vided kgT y|^>o| 2 <? > 200 which could be achieve experi- 
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mentally using Feshbach resonance. 

We acknowledege helpful discussion with S. Nazarenko. 
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